Density fluctuation effects on the exciton condensate in double layer graphene 
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We describe the robustness of an excitonic condensate in double layer graphene against layer 
density fluctuations and the associated charge inhomogeneity, and discuss the implications for ob- 
serving the condensate in current experimental conditions. We solve the mean-field equations for a 
finite imbalance in the Fermi energies in each layer and utilize the results in two phenomenological 
models for inhomogeneity. We find that the stability of the excitonic condensate against density 
fluctuations is strongly dependent on the size of the excitonic gap, and that transport experiments 
(such as Coulomb drag) are promising methods for observing the condensate. 



Excitonic condensates, where pairs of electrons and 
holes in a semiconductor-like system undergo a Bose Ein- 
stein condensation, have been predicted to occur in bi- 
layer semiconductor heterostructuresi They have previ- 
ously been observed experimentally under high magnetic 
fields in the quantum Hall regime* 2 ^ Recent advances in 
fabricating double layer graphene (DLG), where a layer 
of dielectric is sandwiched between two graphene layers, 
have motivated theoretical predictions of the emergence 
of an excitonic condensate at zero magnetic field in these 
heterostructures^ when one layer is electron doped and 
the other is hole doped. In DLG, the thick dielectric pre- 
vents any tunneling of carriers between the two layers, 
while the inter-layer (attractive electron-hole) Coulomb 
interaction drives the formation of an excitonic conden- 
sate. 

The estimate of the critical temperature T c for exci- 
tonic condensation varies over a wide range, even within 
mean-field theory, depending on the level of screening of 
the Coulomb interaction considered in the model. At 
high density, unscreened Coulomb interactions produce 
T c approaching room temperature^ while static screen- 
ing leads to T c < O.IK. 6 ' 7 Large wave vector scattering 
induced by short-ranged disorder further reduces T c j&2 
while the inclusion of dynamic screenin g 10 ' 11 and the full 
band structure^ 2 - tend to favor the pairing. The con- 
troversial nature of the size of the excitonic gap is rein- 
forced by the absence of any signature of the condensate 
in recent Coulomb drag experiments^ on DLG. Since 
disorder scattering strongly suppresses T c j&2 the opti- 
mal regime to search for excitonic condensates in DLG is 
the relatively low-density regime {kpd < 1) where kp is 
the Fermi wave vector and d the inter-layer separation. 
But this effectively low-density (i.e., low fcp) regime in 
graphene is susceptible to strong density fluctuation ef- 
fects (the so-called electron-hole puddles^) arising from 
either extrinsic (Coulomb impurity-induced) or intrinsic 
(e.g., ripple- induced) charge inhomogeneity in the sys- 
tem. In the current work we explore the interesting and 
important question of the effect of inhomogeneous layer 
density fluctuations on the excitonic condensation, em- 
phasizing that this is a distinct physical mechanism from 
impurity disorder,^ where the latter mainly hinders the 
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FIG. 1. Sketch of DLG doped such that the upper layer 
contains electrons and the lower layer contains holes with dif- 
ferent densities. 



condensate formation through momentum-conservation- 
breaking scattering processes. The density fluctuations 
may be caused by any general scalar potential such as 
the Coulomb field of charged impurities or the effect of 
corrugations and ripples. But regardless of the origin of 
the potential, the key distinguishing feature of our work 
is that we allow for imperfect nesting of the Fermi sur- 
faces due to the difference in chemical potential in the 
two layers which is likely to be the primary manifesta- 
tion of disorder in these devices. 

It is well known that the charge landscape of graphene 
mounted on a substrate is inhomogenous^— and that 
when the overall density is small, "puddles" of elec- 
trons and holes are formed. Various mechanisms for 
the formation of this inhomogeneity have been pro- 
posed, including the presence of charged impurities in the 
environmen t 1 7 ' 1 8 and ripples . 19 ' 20 These density fluctua- 
tions can also be viewed as inhomogeneities in the local 
chemical potential iSl The presence of the top and bot- 
tom gates, as well as the dielectric between the graphene 
layers, implies that such chemical potential fluctuations 
are ubiquitous in DLG samples. We study the effects of 
these fluctuations on the excitonic condensate within a 
local density approximation. Keeping in mind the wide 
variation of the excitonic gap in the clean case based 
on the theoretical approximation scheme, we look at the 
two extreme cases: (i) unscreened Coulomb interaction 
and (ii) Coulomb interaction with static screening, with 
a focus on common features obtained by scaling various 
energy scales by Ao, which is the excitonic gap for the 
clean sample in either approximation. 

Our main results are the following: (a) In spatially ho- 
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mogeneous DLG, with different electron and hole Fermi 
energies, the pairing profile depends only on the av- 
erage Fermi energy Jx, and is independent of the differ- 
ence 8[i until \5n\ = 2Ao- Beyond this point, the pairing 
collapses for both the screened and the unscreened inter- 
action, (b) For imbalanced Fermi surfaces, we find no 
evidence of a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
type stated with pairing at finite center-of-mass momen- 
tum |Q| = kp u — kpi in either the unscreened or the 
screened case, (c) We provide two different estimates of 
the effect of charge inhomogeneity on the excitonic con- 
densate: (i) the average pairing gap, which is relevant 
for thermodynamic measurements, and (ii) the fraction 
of the sample area supporting the condensate which is 
relevant for transport measurements. We discuss these 
quantities as a function of the average chemical potential 
and the width of density fluctuations for both screening 
models. Our results show that screening makes the exci- 
tonic condensates more vulnerable to charge fluctuations 
by virtue of the smaller excitonic gap and that transport 
measurements are more robust to these fluctuations than 
bulk thermodynamic measurements. 

The large separation of the two graphene layers implies 
that the local chemical potential (and hence charge) fluc- 
tuations are independent of each other. For convenience, 
we will assume the top layer to be electron doped with lo- 
cal chemical potential fi u and the lower layer to be hole 
doped with local chemical potential — /z;. These quan- 
tities may fluctuate spatially. The independent fluctua- 
tions of the two chemical potentials imply that the perfect 
nesting of Fermi surfaces is lost in these inhomogeneous 
materials even when the global average (fi) is the same 
in each layer. Before going into the details of averaging 
over distributions of chemical potentials through a local 
density approximation, we examine the effect of unequal 
(but spatially homogeneous) chemical potentials in the 
two layers on the excitonic properties. 

The Hamiltonian for a single spin/valley species in 
monolayer graphene in layer A € {it, /} is given by 

H\ v = J2k ( £ k,y - c\-ku c ^ where v € {+,-} de- 
notes the band and — vvk is the single particle en- 
ergy. We will set % = 1 for the rest of this paper. For 
excitonic condensates in DLG, we are interested in the 
case where the average electron (hole) doping in the up- 
per (lower) layer is large (i.e., ^S> A ). Since exciton 
condensation is dominated by processes around the Fermi 
surfaces, we can exclude the filled (empty) valence (con- 
duction) band in the upper (lower) layer from our model. 
This is illustrated in Fig. [TJ The full Hamiltonian is then 
H = H u+ + Hi- + V c h where V e b is the inter-layer inter- 
action between the electrons and holes. We neglect the 
intra-layer interaction since it amounts only to a small 
renormalization of the Fermi velocity in the layerj 23 i 24 
We write the second-quantized operators for electrons in 
the upper layer as etk = c„k+ and for holes in the lower 
layer as 6k = cL_ . This particle-hole transformation al- 
lows us to make a direct connection with the BCS theory 
of superconductivity, albeit with a complicated and nOn- 
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FIG. 2. (a) The pairing amplitude Ak for unscreened inter- 
action with Sfj, = and different p,. (b) The largest pairing 
amplitude Ao for the unscreened interaction at S/j, — 0. (c) 
The evolution of the largest pairing amplitude A with Sfj, for 
the unscreened interaction. The dotted line is Sfi/2 which 
serves as a guide to the eye illustrating the collapse of the 
excitonic pairing when Sfj = 2Ao. (d) Ao for the statically 
screened interaction at Sfj, = 0. 



separable potential kernel. In this notation, the BCS 
pairing ansatz with zero center-of-mass momentum for 
the interaction gives 

kq 

where U(q) is the interaction potential and /(k, k') = 
[1 + cos(#fc — 9k')]/2 is the chirality factor due to the 
graphene band wave functions. Using standard diagonal- 
ization techniques with a mean field in the excitonic chan- 
nel, Ak = (6kfl-k), we obtain two quasiparticle branches 



with spectrum E a ^ — Sfi/2 ± y (vk — fx) + A 2 ., where 
l~l = (ji u + fJi)/2 and Sfi — fi u — /i;. The self-consistent 
gap equation is then given by 



yjivk 1 - /l) 2 + A 2 , 



where the occupation numbers are n a ^{k) = 
<d(—E a (p}(k)) at zero temperature. Equation shows 
that fl and <5/i are the natural variables for analyzing the 
system. We will first look at the system with S/j, = 0, 
where the two Fermi surfaces are perfectly nested. 

Within a model of unscreened Coulomb interaction^ 
which gives the most generous estimate for the excitionic 
gap, the inter-layer potential is V(q) = Vb(q)e~ qd = 
2ire 2 e~ qd /(nq) where V^{q) is the bare two-dimensional 
Coulomb interaction and k is the dielectric constant of 
the environment^ The momentum profile of A^ for 
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Sfi = is shown in Fig. Ufa) for a small inter-layer sep- 
aration of lnm and three different values of p. has 
a non-monotonic momentum dependence with a peak at 
k = kp whose height increases with the size of the Fermi 
surface. This peak height, Ao, is shown in Fig. [2h) as 
a function of p for two different values of the inter-layer 
separation, d = 1, 5nm. The profiles have a convex shape 
which flattens out at large values of p. For Sfi = 0, the 
quasiparticle spectrum E a — — Ep and the positive en- 
ergy branch has a minimum Ao at k — kp. Thus Ao can 
be identified as the true single particle spectral gap in 
the system. 

In the opposite limit, we study the model of Coulomb 
interactions with static screening, which gives the lowest 
estimate of excitonic gaps at high density. In this case, 
the inter-layer potential is given by£ 



oinM; 1 - V b (q) [U u + IL] + V b (q)m u U t (1 - e- 2 ^) 

where H\(q) is the polarization function in a single 
layer— which may be different in the two layers for mis- 
matched Fermi surfaces. The momentum dependence of 
A/j takes a similar form to the unscreened case, with a 
peak at k — kp, but, as shown in Fig.^d), the size of the 
peak Ao is several orders of magnitude smaller. Unlike 
the unscreened case, Ao shows a non-monotonic depen- 
dence on p with a peak occuring at a relatively small 
value of p. The increase in Ao at small p is due to the 
increasing density of states at the Fermi energy, but in- 
creasing the size of Fermi surface also makes the screen- 
ing of the Coulomb interaction more efficient since the 
Thomas Fermi wave vector qtf ~ kp in linearly dispers- 
ing graphene. At large densities, the increased screening 
overwhelms the increasing density of states and the ex- 
citonic gap decreases rapidly.—^ 

We now focus our attention on the effects of Fermi 
surface imbalance on the excitonic properties, using a 
pairing ansatz with zero center of mass momentum. We 
see that changing the sign of Sfi is equivalent to chang- 
ing E a — > —Ep so a discussion of Sfi > will suffice. 
As Sfi increases, we find that the profile of A^ remains 
unchanged until Sfi — 2A . For Sfi > 2A , the pairing 
vanishes for all k leading to the normal state. This be- 
havior is common to both the unscreened and screened 
interaction model and is equivalent to the Clogston- 
Chandrasekhar— limit for loss of pairing in a supercon- 
ductor in a Zeeman field. This is visible explicitly in 
Fig. [U[c), where the peak pairing amplitude A is plot- 
ted as a function of Sfi for the unscreened interaction for 
p = 50 and lOOmeV and d = 1 and 5nm. This behavior 
can be understood from Eq. @ where Sfi appears only 
through the Fermi functions n a rp\. The profile for A& is 
unchanged from the Sfi = case so the occupation prob- 
abilities are unchanged if \5[i\/2 < Ao- Thus the profile 
for Sfi = satisfies the gap equation for Sfi < 2Aq. Be- 
yond this value, the occupation probabilities are changed 
around the gap edge, where the gap edge singularities in 





FIG. 3. Areally averaged excitonic gap, given by Eq. Q. The 
layer separation and screening model are labeled. 



the density of states ensure that the pairing is completely 
lost. 

We have also looked at the possibility of electron-hole 
pairing with finite center of mass in the case of imbal- 
anced Fermi surfaces, specifically for pairing with cen- 
ter of mass momentum \Q\ — \kp u — kpi\, which has a 
rich history in the theory of superconductivity under the 
name "FFLO pairing" .— We have not found any evidence 
of FFLO type excitonic pairing in either the screened or 
the unscreened model. This is not surprising, since, even 
with the model of a constant local interaction, FFLO 
states are known to stabilize in only a narrow parame- 
ter regime in higher than one dimension^ The momen- 
tum dependence of the Coulomb potential would further 
destabilize this fragile state. We safely predict the non- 
existence of any FFLO-typc intcrlayer excitonic superflu- 
idity in density-imbalanced bilayer systems. 

Having looked at the effects of Fermi surface imbalance 
on homogeneous states, we now connect these findings 
to excitonic condensates in actual inhomogeneous DLG 
samples. In the first instance, one could imagine a ther- 
modynamic experimental probe which couples to an area 
of the sample which is larger than the size of the aver- 
age fluctuation (e.g., measurement of a gap derived from 
specific heat). In this case, the measured excitonic gap 
might be best modeled as an average of the gap over the 
spatial area sampled. Since the correlation length of the 
density fluctuations is of the order of the system size, we 
can use a local density approximation to write 

(A) = J J dm dfx u P(fi u )P(m)A(fi u ,fn) (3) 

where P(fJ.\) is the probability distribution for the chem- 
ical potential in layer A. Similar averaging procedures 
have been employed with great success 32 It is well known 
that for monolayer graphene, the distribution of density 
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FIG. 4. Fraction of the area supporting the excitonic con- 
densate in inhomogeneous DLG. The black contour denotes 
A — 0.5 and the white contour is A = 1. 



fluctuations due to charge impurities is Gaussian ; 17 ' 30 
and hence 



2|Ml 



■ exp 



(/M-MIMI) 2 



(4) 



where a n is the width of density fluctuations determined 
by impurity concentrations, (fj) is the global chemical po- 
tential (which we assume to be the same for both layers 
and is set by external gating). It is known from theo- 
retical worki 7 - that in the case of charged impurities, a n 
and the two-dimensional impurity concentration are of 
the same order of magnitude. If the density fluctuations 
are caused by ripples, there is no reason to believe that 
these should be the same in the two layers; on the other 
hand, the field due to charged impurities will be effec- 
tively screened from the more distant layer by the closer 
one. Hence we assume that the inhomogeneities in the 
two layers are uncorrelated. However, inter-layer corre- 
lations can be included by substituting a more sophisti- 
cated function instead of the product P(/i u )P(ij,i). Thus 
we have reduced the effects of inhomogeneity to one phe- 
nomenological parameter a n which takes the same value 
in both layers. Evaluating the integrals in Eq. ^ gives 
the color plots shown in Fig. [3] In the unscreened case 
for small inter-layer distance, the gap is large enough 
that the presence of density fluctuations does not sub- 
stantially alter the gap. However, for static screening, 
the fluctuations are large enough to completely kill the 
exciton condensate for moderate inter-layer separation, 
and even for minimal separation at low density a vanish- 
ing amount of inhomogeneity is required if the gap is to 
persist. 

A more definitive signature of excitonic condensates is 
the superfluid or perfectly inductive response of the sys- 
tem in drag transport measurements. The condensate 



and non-condensate regions form a network, and perco- 
lation theory can be applied to determine the transport 
properties ^ The strongly inhomogeneous samples most 
probably show glassy behavior, but the nature of this 
glassy state requires a more sophisticated treatment of 
disorder. Although a full analysis of percolating clusters 
is beyond the scope of this paper, the fraction of area 
supporting the condensates can be related to percolation 
thresholds. The fractional area of the sample which sup- 
ports the condensate is given by 
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dn u dfiiP(fx u )P(fXi)e(fi utM )e(2A - \Sfx\) (5) 



The first step function represents the fact that in the in- 
homogeneous case it is possible for /i u and /i; to be in 
the same band. If this is the case then the condensate is 
not allowed. Figure @] shows A for both the unscreened 
and statically screened interactions, and for the two val- 
ues of the inter-layer separation shown previously. The 
black contour denotes A = 0.5, where half the area sup- 
ports the condensate, and the white contour shows A = 1 
meaning that the whole sample is excitonic. We see that 
in the unscreened case, the condensate should be sta- 
ble at experimentally accessible values of density fluctu- 
ations, but in the statically screened case the fluctuations 
have to be smaller even than those currently available in 
hexagonal boron nitride (hBN) devices! 30 ' 32 By compar- 
ing Figs. @] and El we see that in samples with density 
fluctuations, the formation of the excitonic condensate 
is easier to detect in transport-style experiments than in 
probes that average over the bulk of the sampled 

We also note that other author o 10 ' 11 have considered 
the effect of dynamic screening on the exciton conden- 
sate without disorder. In the high-density limit (where 
static screening gives a tiny gap due to the large size 
of the polarization function), dynamic screening gives a 
gap which is closer to the unscreened interaction than the 
static screening. Using the data presented in Ref.[I3, we 
estimate that for physically realistic parameters, the dy- 
namic screening theory predicts that the peak gap in the 
high-density regime will be ~ O.lmcV, indicating that 
the numerical effect of disorder is roughly similar to that 
of the static screening in the low density regime. Since 
our analysis shows that the nature of the screening has 
no impact on the effects of inhomogeneity, Figs. [3] and 
2] show that a gap of this size is marginally observable. 
Therefore, just as the low density statically screened in- 
teraction may be observed in the cleanest systems, so 
should a gap generated by dynamic screening in the high- 
density regime. 

In summary, we have examined the impact of chemical 
potential fluctuations and charge inhomogeneity on the 
formation of an excitonic condensate in DLG. The origin 
of these charge fluctuations may be from ripples in the 
graphene surface, charged impurities in the environment 
of the graphene, or any other general spatially-varying 
scalar field. While the absolute value of the gap depends 
sensitively on the model of interaction used, the effects 
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of Fermi surface imbalance are identical in the two cases. 
However, since the inhomogeneity can be controlled only 
on an absolute scale, we find that the size of the gap is 
the crucial factor which determines the robustness of the 
condensate against this form of disorder. If the gap is 
of the order of a few meV, then the level of fluctuations 
found with hBN substrates should be sufficient to observe 
the condensate in Coulomb drag experiments. Finally, if 
static screening is important in the formation of excitonic 
condensates, there is an optimal Fermi energy (which is 
~ 20meV for d = lnm) , where it would be easiest to see 



the condensate. 
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